Note on the ring approximation in nuclear matter 

E. Bauer* 

Departamento de Fisica, Universidad Nacional de La Plata and 
Instituto de Fisica La Plata, CONICET, 
C. C. 67, 1900 La Plata, Argentina 
■(^" ■ (Dated: May 1, 2008) 

o ■ 

S ! Abstract 

The response function to an external prove is evaluated using the ring approximation in nuclear 
matter. Contrary to what it is usually assumed, it is shown that the summation of the ring series 
and the solution of the Dyson's equation are two different approaches. The numerical results 
exhibit a perceptible difference between both approximations. 
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The ring approximation is widely used in many nuclear physics problems. It consists 
of the infinite sum of one particle-one hole bubbles, where the Pauli exchange contribution 
is neglected [1]. In fact, the ring approximation is the direct part of the Random Phase 
Approximation (RPA). By neglecting the exchange terms, the particle-hole series reduces 
itself to a geometric series, which is easily summed up. As an alternative derivation, it is 
usually proposed the ring series as the solution of the Dyson's equation [2] (see also [1]). It 
has come to us as a surprise the existence of an inconsistency between the interpretation 
of the ring approximation as a solution of the Dyson's equation and the explicit evaluation 
(employing the Goldstone's rules) of the particle-hole diagrams which originates the ring 
series. This inconsistency comes into play only when one particlc-holc configuration is on 
the mass shell, for example, when we study the response function. In this contribution we 
discuss this point and we show that there are two alternative approximations. One is the 
solution of the plain Dyson's equation and the other is the summation of the ring diagrams. 

Let us consider an arbitrary one-body operator 0^ exciting the nucleus from its ground 
state |0 >. The action of Ca is characterized by the nuclear response function per nuclear 
volume, 

'5a(go,q) = -\ Im < Q\OjG{qo)Oa\0 >, (1) 
where go and q are the energy and momentum transfer, respectively and Q is the nuclear 
volume. The nucleon propagator G{qo) is given by, 

1 1 
^ qo-H + ir, - qo + H + ir, 
being H the nuclear Hamiltonian. As usual H — T + V, where T is the kinetic energy and 
V is the residual interaction. The identity is expressed as, 

/ = ^|n><n|, (3) 

n 

where \n > represent a compleat set of orthonormal states of H. By inserting the identity 
twice in Eq. (1), we have, 

Sa{qo,q)^ ^Yl \<n\Oa\0>\^5{qo-qn), (4) 



n 



where qn = E^ — Eq {h — c— 1), and E^ are the excitation energies of the eigenstates \n >. 
In the present contribution we explore three different external proofs, 

0^ = Y,e^^-^^d^^,) (5) 



with Xj denoting the intrinsic coordinate for individual nucleons and the sum j runs over 
all nucleons. In this equation, — 1; ^L{j) — cr{j) ■ Q'^z{j) ^-nd Cr(j) = cr(j) x Q'^z{j) 

which are usually named as the isoscalar central and isovector spin-longitudinal and spin- 
transversal operators, respectively. The next step is to propose a model for the Hamiltonian. 
The simplest choice is to keep only the kinetic energy T. In this case, the final state is a one 
particle-one hole excitation as the one drawn in the first diagram in Fig. 1. Before we show 




FIG. 1: Goldstone's diagrams representing the firsts contributions to the ring approximation. In 
each diagram an up (down) arrow constitutes a particle (hole) , a wavy Une is the residual interaction 
and a dot stand for the external operator. It has been added an horizontal dashed line to indicate 
the configuration on the mass shell. 

the response function, in Appendix A, it is presented the lowest-order polarization insertion 
n°(go, q), which is further expressed as the sum of it real and imaginary parts, 

U'^{qo,q)=n{qo,q)+tI{qo,q) (6) 
Now the response function to Oc, is, 

SUqo,q)^I (7) 

Next, the residual interaction V is incorporated. As a model for V, we use the one given 
in Eq. (A2), which can be rewritten as a sum of a isoscalar central and isovector spin- 
longitudinal and spin-transversal terms (see Eq. (A4)). One way of taking care of the 
residual interaction is by means of the Dyson's equation, where a higher-order polarization 
insertion is obtained as, 

W'y^qo, q) = n°(go, q) + n°(go, q) V{q) n^^^(go, q). (8) 
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In this equation the Pauh exchange terms have been aheady neglected. For this reason, this 
is an algebraic equation which solution is the sum of a geometric series in VU^, 

Using the solution of the Dyson's equation a new response function is obtained. This is done 
by replacing X in Eq. (7) by /mll^^*, 

sSHtM Q) - (l-Vc(.,T)7^)^ + (Vc(.,T)X)^■ 

The Eq. (9) is usually interpreted as the sum of a series of one particle-one hole bubbles. 
The firsts terms to this series are shown in Fig. 1. In these diagrams, a horizontal dashed-line 
indicates that this configuration is on the mass shell. It is interesting to analyze each term 
in the series separately, we expand Eq. (9), 

^-i^^ = n°(i + vu' + {vur + ...), (11) 

by taking the imaginary part of this sum, we have, 

/m(n°) = I 
7m(nVn°) ^ 2nvi 
7m(n°(yn°)2) = (37^'2: - 2:=^) 

(12) 

By inspection of the diagrams in Fig. 1, each term has the following interpretation. The 
first term (zeroth power in V), is the first diagram in the left hand side in this figure. Using 
the Goldstone's rules, the analytical expression for a one particle-one hole bubble is given 
by n°. When the bubble is put on it mass shell, we take the imaginary part, X. The 
next contribution (first power in V), is represented by the second and third diagrams in 
the same figure. In the second (third) diagram the lower (upper) bubble is on the mass 
shell. Analytically, the bubble on it mass shell is given by X, while the other bubble (in the 
same diagram) is off the mass shell, TZ. As both contributions (second and third diagrams) 
are identical, one has a factor two (i.e., 2TZ V I). This association between diagrams and 
physical states fails for the next order contribution. There are three contributions, where 
the first one is shown in Fig. 1, while the two remainders ones are the same draw, but 
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with the dashed-hne (which represents the configuration on it mass shell), in the middle 
and upper bubble. The first term (i.e., ?>R?X) is easily interpreted as the sum of these three 
contributions, where only one bubble at a time is on the mass shell and a factor three results 
from the equality of the three contributions. However, the X^-term can not be interpreted: 
in terms of Eq. (4), all diagrams represent the square of a transition amplitude. To put it in 
other words, in one diagram only one configuration can be on the mass shell. The X^-term 
would imply a diagram with three bubbles simultaneously on the mass shell. 
We go back to Eq. (9) where we keep only the terms compatible with Eq. (4), 

7m(n°) = X 

7m(nVn°) = 27^yx 
/m(n°(yn°)2) = 37^^ v^x 

7m(n°(yn°)^) = (Ar + i)7e^y^T 

(13) 

Each term in this series has a straightforward physical interpretation in terms of the so- 
called ring diagrams. For this reason, we call the sum as H''*"^. The summation can be 
easily performed once we notice that, 

d ( 1 



d{nv) \. 



1-TZV 
The final result for the sum is. 



1 + 2nv + 3{nvy + A{nvy + . . . (u) 



It should be noted that the expression in the left hand side in Eq. (11) is the sum of the 
series in the right hand side, as long as |y n°| < 1, therefore we have. 



|yn°| = \V\VWTX^< l, asn,XeRe^ \Vn\ < l.-. l-Vc(L,r)^7^0. (16) 

Therefore, if the sum in Eq. (11) exists, so does the one in Eq. (15). In order to obtained 
gDys g^j^^ grmg^ ^^^y the imaginary part in Eq. (9) has been evaluated. For completeness, 
in Appendix A the real parts are also calculated. 
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As mentioned, we have two different response functions, S^''^^ and S"'^^. In Fig. 2 we have 
plotted the numerical result for these two functions, just to show that the difference between 
them is relevant. It should be emphasized that both response functions are valid solutions 
for two different ways of dealing with the residual interaction V: the S'^^*-response is the 
solution of the Dyson's equation and the S'""^-response is the sum of the ring diagrams. In 
the present contribution, it is claimed that the interpretation of the solution of the Dyson's 
equation in terms of ring diagrams is wrong, as long as the polarization insertion 11° has a 
not-null imaginary part. In some physical problems, such as the study of zero sound [3, 4] 
or core polarization [5], only the real part in Eq. (11) is needed. By inspection of Eqs. (A5) 
and (A8), it is easy to check that U^^' = W'"'^ when 1 = 0. In this case, the solution 
of the Dyson's equation is also the sum of the ring series. This element could have been 
misleading in the former interpretation of the solution of the Dyson's equation. Before we 
end this paragraphs, another point should be addressed: the interpretation of the solution 
of the Dyson's equation in term of Eq. (4). This can be done as follows. We work with the 
residual interaction, 

V^'^qo, q) = V{q) + V{q) n°(go, q) T^^^^go, q)- (17) 
which has the simple solution, 

yDys ^ K 

When this complex interaction is used in replacement of V in the second and third dia- 
grams in Fig. 1, a solution of the Dyson's equation compatible with Eq. (4) is obtained. 
Analytically, using Eqs. (9) and (18), it is straightforward to check that, 

n^2^* = n° + nV^^^n". (19) 

In this case, only three Goldstone's diagrams comes into play (the first, second and third 
in Fig. 1, with the physical states as marked in this figure). An expansion in term of ring 
diagrams of Eq. (18) is possible, but making no connection with physical states. 

As a further quotation, the Dyson's equation can be split into it the real and the imaginary 
part. In any case, the solution is the one given in Eq. (A5). If II"'*^ (see Eq. (A8)), is replaced 
in the Dyson's equation, the imaginary part of this equation is satisfied, but not the real 
part. This observation is given only as a warning: in the cases in which the Dyson's equation 
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FIG. 2: Response function per nuclear volume. In each graph, the dot line represents S^, the 
continuous one S"'^^ and the dash-dot S^^^. The momentum transfer by the external operator 
is chosen as q = 400 MeV/c, while the parameters entering in V are / = 0.3, g' = 0.5, A^r = 
1300 MeV/c and Ap = 1700 MeV/c. For the Fermi momentum it has been used, kp = 1.36 1/fm. 
The functions Sc,lt are given in units of MeV^i fm'^ xlO"^. 

is solved numerically, no matters if it is needed only the imaginary part of the solution. Both 
real and imaginary parts should be found. 

As a concluding remark for this contribution, we have discussed the response function 
employing two different ways of dealing with the residual interaction. The first one is by using 
the solution of the Dyson's equation and in the second, we have analyzed the ring diagrams. 
For the ring diagrams, we have taken special care of the configuration which is on the mass 



7 



shell and the interpretation of these diagrams in terms of the Eq. (4). A similar analysis for 
the solution of the Dyson's equation has been proposed for completeness. Both analytically 
and numerically, these solutions are different. They represent different approximations and 

they are both correct. A step forward in this kind of analysis would be the discussion of 
the Continued Fraction Approximation [6], a subject which has been paid some attention 
recently [7]. 
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APPENDIX A 

The lowest-order polarization insertion is, 

n°(5o, q) = 4 f^,e{\q + k\- k,)e{kF - k) ( — ^ ^ ) 

J K^T^Y \qo - tq+k + tk + IV qo + tq+k-tk-tvj 

(Al) 

where tp — p^/(2m), with m being the nucleon mass. In this equation kp is the Fermi 
momentum. 

We present now our model for the residual interaction V, 
V{q) = ilrl{q){f + g'(T-(T'T-T' + V^{q)a q(T' qT T' + V,{q){(rxq)-{(T'xq)T-T'), (A2) 

A*7r 

where it has been taken the static limit. Therefore, V.,r{q) = —(f/io.^ + /^^) and Vp{q) = 
— (Tp/TT^YCpq^/{q^ + /x^), where /i^ (/Xp ) is the pion (rho) rest mass, fl/^-n = 0.081 and 
Cp = 2.18. The form factor of the ttA^A^ {pNN) vertex is (F^), where F^ = ((A^ - 
m^)/(A| -\-q^)Y. Using the property, 

a ■ a' ^ a qa' q + {a x q) ■ {(t' x q), (A3) 

the Eq. (A2) can be rewritten as. 



V{q) - 4^l{q){Vc + Vl(T ■ qa' qr -r' + Vt(o- x q) ■ {(t' x q) r ■ r'), (A4) 

with obvious definitions for Vc,l,t- 

As a final point for this Appendix, we split the solution of the Dyson's equation (Eq. (8)), 
into it real and imaginary parts, 

^ n{i-vn)-vi^ I 
{i-vnf + {viy {i-vny^ + {vTY ^ ' 
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We now perform the same procedure as in Eq. (13), but for the real part of the ring series, 

i?e(n°) = n 

i?e(nVn°) ^ n^v 

i?e(n°(vn°)2) = n"" 



where the sum is. 



Finally, we can write. 



N 



(A6) 



i?e(n""») = — - — . (A7) 



- {1-vny + {i-vnr' ^^^^ 
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